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1 Introduction 

Machine issues at particle factories are dominated by luminosity optimiza- 
tion, which is the overlap of the density functions p of the two beams over 
space and time. For a single beam crossing, 

L = c J dVdTpi(r,T)p 2 (r,r), (1) 

where dV is the volume element and r is a time variable of order of the beam- 
crossing duration. Optimal luminosity is achieved by perfect transverse 
overlap of two equal and Gaussian beams squeezed to the limit allowed by 
the tune shift. For a single beam crossing, that reads 

Lo{t) = ^ N2{ f\ C , r I ^dre-^V^W+J/V^W+C^+C-) 2 )/-^*)) 
uv ; (2-afa x {t)a y {t)(x z {t) 
JVi(t)JV 2 (t) 



where the N\ % and &x,y,z are the beam populations and spatial dimensions 
at any given time t (which is a run-time variable of order one hour) . 

This formula becomes rather complex over time. Particles are deflected 
by the other beam at each crossing, significantly affecting the Twiss matrix 
of the machine. The beam currents iVj(i) decrease due to beam lifetime also 
causing the machine's Twiss matrix to drift. If the machine is perfectly sym- 
metric, the transverse dimensions will change but the beams will maintain 
perfect overlap. 
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Even symmetric machines have some degree of asymmetry, and beams 
start moving independently in the transverse plane as soon as the run starts. 
At .B-Factories such as CESR, PEP-II, and KEK, beams have horizontal 
dimensions a y of order several microns, with aspect ratios e = cr y /cr x ~ 
0.02 — 0.04. A drift of 5 microns is enough to spoil the luminosity. 

A better description of the luminosity over time should be 



w(t) is the positive-defined waste parameter due to non-instantaneous opti- 
mal overlap. If w(t) is known, the wasted integrated luminosity is defined 
as 



where / is the machine frequency. The waste parameter can be readily 



The waste parameter is clearly of great interest, and one of the most im- 
portant issues not only at B factories but also at future linear colliders. As 
soon as w is non-zero, a correction should be applied to restore optimal lu- 
minosity. The ability to measure and reduce w would make for a substantial 
increase in the delivered luminosity of any machine. 

Although w can be defined mathematically from Equation the beam- 
beam topology cannot be measured directly. Techniques have been devel- 
oped that measure the transverse displacement of the centers of gravity of 
the beams by BambadeQ and Sagan, Sikora and Henderson Q. Both of 
these techniques actively displace one of the beams, and monitor the other 
beam to observe the strength of the beam-beam interaction. In practice 
these techniques are sensitive to the relative displacement of the beams cen- 
ters from optimal beam-beam overlap. 

Generally, a discussion of the waste parameter must include all possible 
degrees of freedom in the evolution of a machine over a run. There are 
seven parameters which can affect optimal beam-beam overlap Q. These 
are shown in Figure |]. Briefly there is a transverse displacement between 
the two beam centers described by a vector (Ax, Ay), the transverse sizes 
of the two beams a x i, a y %, a X 2, and a y 2, and a relative rotation of the two 
beams <f>. The two beam currents also affect the beam-beam interaction. 
Currents are easily monitored and are not included in the discussion below. 



L(t)=L (t)(l-w(t)). 



(3) 




(4) 
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In this paper, a technique is proposed by which six of the seven pa- 
rameters can be passively monitored with the observation of wide angle 
beamstrahlung. In case of non-zero waste, which is called a "pathology", 
the responsible parameter is identified unambiguously, and the amount of 
needed correction is measured. The seventh parameter can easily be mea- 
sured in a beam scan also using the wide angle beamstrahlung signal. 

Seven parameters to characterize the beam-beam collision is a large num- 
ber. It is easiest to discuss the problem if it is broken into two parts. 

• The machine is perfectly symmetric, that is, the machine optics is ex- 
actly the same for both beams. In Figure [l], that means that the two 
beams have zero offsets, zero rotation, and the same transverse dimen- 
sions, resulting in only two parameters. Dropping indices, they are the 
transverse dimensions a x and a y . If the machine is symmetric, beams 
maintain optimal overlap, but the optics is affected by the varying 
currents. The luminosity is determined by the transverse size of the 
beam. The case of measuring the transverse beam size is discussed in 
Section 3. 

• The beams move independently in the transverse plane due to ma- 
chine asymmetry decreasing overlap and luminosity. In section 4 the 
measurements of the relative sizes of the two beams, their transverse 
displacement, and the angle between them are described. 

In this paper the use of large angle beamstrahlung, which is described in 
detail in reference^], is described as a beam-beam monitor that allows com- 
plete control over both the beam-beam interaction strength and transverse 
displacements. Large-angle beamstrahlung observables, combined in a sim- 
ple 2-dimensional diagram which is called the beamstrahlung diagram, mon- 
itors the wasted luminosity. 

In Section 2 the information content of large angle beamstrahlung is 
discussed. Section 3 covers the symmetric machine case, concentrating on 
measurements of the beam size. Section 4 covers asymmetric machines, in- 
troduces the beamstrahlung diagram, and shows how the waste parameter 
can be measured. The use of the beamstrahlung diagram to eliminate wasted 
luminosity is shown in Section 5. Three appendices are included for com- 
pleteness. Appendix A derives in a simple way three crucial properties of 
large angle synchrotron radiation. Appendix B provides a description of the 
beam-beam simulation developed for this paper and Appendix C evaluates 
the simulation's accuracy. 
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2 Large Angle Beamstrahlung 

The properties of large angle radiation, emitted by a ultra-relativistic parti- 
cle, differ dramatically from the classical synchrotron radiation formulae [0]. 
Appendix A shows that the approximations used in reference || and in this 
paper are valid at large angles for all present and proposed e + e~ collid- 
ers, if beamstrahlung detection is to be done at or near the "magic angle" 
described in Q. Three properties of large angle radiation are derived in 
Appendix A. Of particular interest is the 100% linear polarization, either 
parallel or perpendicular to the bending force, obtained at certain azimuthal 
locations at large angle. 

At CESR for example, it is possible to detect such radiation in vis- 
ible light at a location 5 meters away from the interaction point, at a 
6mrad angle. The beam-beam interaction occurs over a volume of order 
300/iinx7/mix7mm, and particles are typically deflected laterally by 10~ 2 
mrad. Thus the light detector is seen at the same angle by all of the beam, 
and throughout the dynamic beam-beam collision. These are the conditions 
termed as "CESR conditions" and used for the calculations of Sections 4 
and 5. A fixed fraction of the beamstrahlung energy is collected at such 
a location, effectively measuring the total energy up to a constant. Dif- 
ferent polarization components can also be easily observed, by filtering the 
observed light through polarimeters. 

The two polarization components can be used to build the radiation 
vectors Ui from one beam and U2 from the other beam, which are two- 
dimensional vectors in the first quadrant. The first dimension is the horizon- 
tal component of the polarized beamstrahlung power signal and the second 
is the vertical. The total energy vector U is defined as Ui + U2. At large 
angles the polarization components and radiation spectrum factorize|J and 
a different orientation of the polarimeters would simply rotate the horizontal 
and vertical axes.. 

As mentioned in the Introduction, at present and proposed machines, 
beams are very flat (e ~ 0.02 — 0.04). It is convenient to develop the theory 
only for flat beams which leads to two simplifications. First, terms of order 
e and higher can be neglected in equations as needed. Second, a natural 
preferred orientation exists in the transverse plane, which is adopted to 
produce the results of this paper. 

It should be noted that two counters on each side, each looking at a dif- 
ferent polarization component, and in absence of background, are enough to 
extract complete information from beamstrahlung. As an example, given the 
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formulae in Appendix A, U x can be measured by measuring the x— polarized 
component at zero degrees in azimuth, and U y by the x— polarized compo- 
nent at 45 degrees. 



3 Symmetric Machines 

If a machine is perfectly symmetric, the beam currents and transverse di- 
mensions of the beams will evolve, while maintaining perfect overlap. Mea- 
surements of the beam sizes a x and a y , determine the luminosity. In this 
case most of the interplay between machine and beam-beam interaction is 
through the dynamic beta effect. 

The dynamic beta effect is proportional to the average electric field seen 
by one particle over many beam crossings, hence it is proportional to the 
charge in the other beam, times the average inverse impact parameter b 
between particles of beam 1 and particles of beam 2|jJ, 

< Ei >oc N 2 < -r > . (6) 
fr 

From|| the beamstrahlung energy is proportional to 

C/i oc Ni < E\ >, (7) 

or 

< E\ >oc Ux/Ni. (8) 

The < Ei > and < E\ > are related through the transverse shape of 
the beam, which can be taken to be Gaussian with no loss of precision. 
Therefore, monitoring the dynamic beta effect can be done efficiently by 
monitoring the U{ and the Ni at the same time. 
Equation ^ can be rewritten as[||, ||] 

Ui oc ^f(e). (9) 

The beam length, a z , is usually constant, and will not be considered here, 
but clearly a beamstrahlung detector can also be used to monitor the beam 
length, for example during machine studies. The function /(e) varies slowly 

/(e) ~ 1 + 11.4c, (10) 



and can be considered nearly constant in the following. 
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The result above assumes "stiff" beams. A stiff beam is one where the 
beam particles do not change their transverse position appreciably during 
the collision. Appendix B shows that dynamic effects are negligible. 

In flat beams most of the impact parameter is due to the distance in x be- 
tween the particles, and the energy radiated is almost only dependent on a x . 
For perfect overlap of stiff Gaussian beams the energy U is unpolarized||. 
No information can be extracted out of polarization, and beamstrahlung 
cannot monitor passively symmetric changes in a y . The total power radi- 
ated is thus sensitive to a x . 

However, as pointed out in references [||, |(|, a scan of one beam along the 
vertical axis will produce the characteristic camelback feature in the plot of 
U versus the beam-beam offset seen in Figure ||, which has already been used 
in the detection of beamstrahlungf8). The o y can be precisely determined 
by measuring the peak-valley distance d shown in Figure ||. The relation 
between d and o y is 

d ~ 3.97cr y (l - 5.4e). (11) 

Currently, the CESR beams are artificially perturbed with an amplitude 
of order 0.01a y to measure the beam-beam interaction by observing the 
effect of the perturbation on the other beam via the lock-in effect^]. It is 
conceivable that this technique could ultimately be used to determine a y 
without scanning. Note that a beam scan could also be used to measure a x 
separating it from a z . 

A beamstrahlung monitor can be very useful even when a machine is 
perfectly symmetric, allowing purely passive monitoring of the beam-beam 
interaction and thus the beam length or a x . It can be used to measure 
a y in a beam scan. The next Section, which deals with purely asymmetric 
pathologies, shows that this method is truly valuable when beams are not 
colliding head on in the transverse plane and may have different transverse 
sizes. 

4 Asymmetric Machines 

If a machine is asymmetric, as all real machines are to some degree, the 
two beams will drift independently in the 7-dimensional space that induces 
luminosity waste. For the purpose of studying asymmetric machines, a single 
pass beam-beam simulation program was written. The program generates 
complex beam-beam configurations involving all the pathologies shown in 
Figure |l[ These configurations are, in principle, computable analytically in 
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Paramter 


Symbol 


Value 


Beam Width 




300 /jm 


Beam Height 


Oy 


7 fim 


Beam Length 


<Jz 


1.8 cm 


Bunch Charge 


N 


11 x 10 10 


Relativistic Factor 


7 


10 4 



Table 1: Beam parameters chosen for the simulation results presented here. 



the limit of stiff beams. It was important also to cross-check the effects of 
beam-beam dynamics, as the particles of one beam are deflected towards 
the center of the other beam. The latter is an effect that must be computed 
by simulation. 

The simulation program is described in Appendix B. Its precision is eval- 
uated in Appendix C and is found to be between 0.1% and 0.2%, for beam- 
strahlung yields, and better than 1%, for the luminosity enhancement due 
to beam-beam pinching. The nominal conditions to produce results in this 
paper were chosen as in Table |] and are appropriate for l?-Factories. The 
measurement of the two beam sizes was discussed in the previous section. 
The remaining five parameters are discussed here: two relative transverse 
dimensions; two offsets; and one angle. Present day beam position monitors 
have spatial resolutions of order 20 /jm, which is substantially less than the 
a x of these beams, and should always provide adequate overlap along the 
x— axis leaving four pathologies of concern. An offset in x will generate a 
unique configuration of the beamstrahlung diagram, which mirrors the one 
obtained for an offset in y which is discussed below, and can be analyzed in 
a completely equivalent way. 

For simplicity it is assumed that only one beam is developing a pathology 
at any given time. The four pathologies that lead to wasted luminosity 
are shown in Figure [|. They correspond respectively to a vertical offset, 
imperfect vertical focusing, imperfect horizontal focusing, and a rotation of 
one beam with respect to the other. All these pathologies except the third 
have been observed at CESR. 

The expected value of each polarization component is also assumed, for 
optimal beam-beam collision, which we call Uq. In practice Uq can be mea- 
sured experimentally by continuous observation of colliding beams, or by 
separately determining the beam currents, and a x and <r y with a beam-beam 
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scan as discussed in the previous section. 

The beamstrahlung diagram plots Ui,U2 normalized by Uq. In the 
figures below the contribution from the pathological beam is represented 
by the dashed arrow. The diagram has four degrees of freedom. The total 
power monitors the beam-beam interaction strength, and three independent 
dimensionless asymmetries can be defined. 

As mentioned in Section 3 if the collision is perfect and the beams are stiff 
the beamstrahlung radiation is unpolarized. Thus the normalized U's are 
equal and the vector from each beam in a perfect head-on collision are on top 
each other at 45 degrees as shown in Figure ^. With the Uq normalization 
one obtains the perfect collision point at (1, 1) for both beams. 

The effect of dynamic beams can be estimated from Table |2| in Ap- 
pendix C. For example at CESR dynamic beams increase U x by 0.9% and 
U y by 2.7%, moving the perfect collision axis 0.5 degrees above 45 degrees. 
Such a small modification is nearly invisible in Figure ||| and can be neglected. 

Figure 5 shows for stiff beams the beamstrahlung diagrams for each 
pathology shown in Figure 3. Each has a unique pattern, which a feedback 
algorithm can discern and correct. In general, if beam 1 is smaller in x(y) 
than beam 2, then it will radiate less energy in x(y). 

Figure 6 is the same as Figure 5, but for dynamic beams. Comparison 
of the two figures shows very little difference. The effect of dynamic beams 
is small. Thus the beamstrahlung diagram presented in this paper is a 
universal display of the pattern associated with beam-beam pathologies at 
CESR, PEP-II, KEK, and in the future at a ~1 TeV e + e~ machine. 

Figure 7 is the same as Figures 5 and 6, but with an offset in x, 18/rni, 
or 0.06(72, comparable to the resolution of beam position monitors. Again 
very little change is observed with respect to Figure 5 showing that small 
horizontal offsets have small impact. 

Asymmetries corresponding to each pathology in Fig. 3 are defined as 



where O is the Heaviside function meaning in this case that the asymmetries 
A{ are not defined when the argument of the Heaviside function becomes neg- 
ative. The indexing was chosen to indicate that the second, a beam bloated 
vertically, and third, a beam bloated horizontally, pathologies are generated 




(12) 
(13) 
(14) 
(15) 
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from both a zero dipole moment and a non-zero quadrupole moment in the 
transverse charge distribution, and as such they should be equally ranked. 

These asymmetries are not independent. The usefulness of these beam- 
strahlung asymmetries is shown in Figure 8 which displays their dependence 
on the waste parameter defined in Equation [3|. Each asymmetry's contribu- 
tion to the waste parameter of Section 1 is 

dw 

- ~ 8A A - < 16 > 

but they can not be summed together because they are not independent. 

Evidence is provided in the next Section that the asymmetries should 
be minimized strictly in the order defined by Eqs. 12-15. In short, the total 
waste parameter can be defined as 

w ~ ^2 dwiAt, (17) 

i 

where the hatted derivative is defined as 

£>Wi = {^-)A j =min.,j<i- (18) 

Eq. 17 represents the main result of this paper. The derivatives are com- 
puted, the asymmetries are measured, and the waste parameter is obtained. 
Note that if the asymmetries were completely independent, the specifications 
Aj = min. would have been unneeded. Asymmetries 2 and 2' represent both 
quadrupole corrections, and can be interchanged without harm. 
For horizontal offsets between the two beams an asymmetry 

At = (U X /Uy - l)e(U X /Uy - 1) (19) 

can be defined. 

We note that for a 10% change in luminosity, the values of the asymme- 
tries change by 0.1 for A\ and A^, 0.25 for A<i and A' 2 and 0.05 for ^3. Thus 
these asymmetries have excellent sensitivity to wasted luminosity. 

5 The Virtual Operator 

Here examples are shown of how the beamstrahlung diagram and the asym- 
metries defined in Equations 12-15 and [l^ can by used to eliminate wasted 
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luminosity even in the presence of multiple pathologies in the beam-beam 
collision. 

We demonstrates this by studying the complete set of six double patholo- 
gies, shown in Figure 9, which can be derived from the four single pathologies 
shown in Figure 3. Figure 10 represents the beamstrahlung diagrams cor- 
responding to the pathologies displayed in Figure 9. A feedback program, 
dubbed the Virtual Operator, finds the highest-ranking asymmetry, mini- 
mizes it by changing the appropriate collision parameter, and obtains the 
beamstrahlung diagrams of Figure 11, displaying only one pathology which 
is trivial to correct. 

Two comments are in order. First, if the largest asymmetry were to 
be corrected first, instead of the highest-ranked, convergence would not be 
achieved. To prove the point, it is enough to compare Figs. 5b, lOd and lid. 
If the sextupole correction is acted upon first, A3 would have to increase as 
opposed to being minimized. 

Second, we wish to prove that minimization of a higher-ranked asym- 
metry effectively corrects the associated pathology. Although all the double 
pathologies were tried, only Fig. 9c, which corresponds to a vertical offset 
plus a rotation, is presented. Ai and ^3 are the two most correlated asym- 
metries. The asymmetry A\ is not zeroed, and cannot be zeroed by moving 
one beam. Figure 12 shows the dependence of A\ versus the vertical offset, 
showing that minimization of the asymetry gives the desired correction. 

We did not consider horizontal offsets here, but they could easily have 
been included, as discussed in the previous Section. 

6 Conclusion 

The beamstrahlung diagram and asymmetries derived here demonstrate a 
complete and rigorous method for luminosity optimization. The wasted 
luminosity is for the first time related to quantities that are instantaneously 
observable, and specify the necessary correction. We have considered a 
complete class of beam-beam pathologies. 

If the machine is perfectly symmetric a beamstrahlung monitor is very 
useful for measuring the size of the beam. In the case of asymmetric beams 
a beamstrahlung monitor is extremely powerful. The study of the beam- 
strahlung diagram derived from the power and polarization of the beam- 
strahlung signal allows identification of the beam-beam pathology, identifi- 
cation of the "bad" beam, and measures the correction that needs to be ap- 
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plied. In short the wide angle beamstrahlung signal analyzed in the manner 
described here is a powerful tool to eliminate wasted luminosity at present 
and future e + e~ colliders. 

Appendix A 

The properties of short magnet radiation were first discussed by Coisson in 
reference ||]. In the classical model Q, the bent electron is made to sweep 
through the detector in a "searchlight" fashion, effectively covering all beam- 
detector angles. In the Coisson's model the opposite extreme is adopted, and 
the angle is kept constant throughout the orbit, the large angle approxima- 
tion. Both models predict the same power, the same total polarization, and 
the same typical angle, of order I/7 for the emitted radiation, but they differ 
dramatically in the spectrum at large angles. 

The Coisson model is of interest here because the detector's angle is con- 
stant throughout the collision at colliders such as CESR. At large angle the 
classical model predicts a steep fall-off of the power, exponential both in the 
photon energy and in the cube of the observation angle. The Coisson model 
predicts three properties of large angle beamstrahlung radiation. They are: 

• The cutoff energy, at large angle, does not depend on 7. There is 
no exponential fall-off as predicted by the "searchlight" approxima- 
tion, making detection possible. In particular at 6 mrad at CESR, for 
example, visible radiation is at or below the cutoff frequency. 

• The polarization is linear at a fixed location in azimuth with an eight- 
fold pattern, (cos 2 2(f), sin 2 2<^>) around the azimuth. The angle <f> is the 
angle between the net transverse force experienced by the beam and 
the detector location. Thus the pattern of the polarization provides 
information about the beam-beam overlap. 

• The large angle double differential spectrum is proportional to (7$) , 
and not exponential. The large angle power scales as I/7 2 . Thus 
the situation at B factories is more favorable than at higher energy 
machines. 

These properties are re-derived here in an elementary way for constant 
large angle of detection. Consider an extremely relativistic particle, 7 >> 
1, undergoing a vertical deflection, due for example to a horizontal dipole 
magnet exerting a force F over a length a z . Radiation of energy k = huj 
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is detected at an angle 9 which is much larger than 1/7. In the laboratory 
frame the radiated energy is equal to § 

U=\^F*a z . (20) 
3 mc A 

A simpler derivation is possible by studying the radiation in the rest frame 
of the radiating particle. Note that all quantities in the particle rest frame 
are starred as shown in Figure 13. The radiation will have a dipole pattern 
with angular intensity proportional to the squared sine of the angle between 
the direction of detection and the direction of the force. The force maintains 
its vertical direction and has a modulus 

F* = 2-fF. (21) 

The angle is very large in the laboratory frame, and the corresponding di- 
rection in the rest frame is very close to the backward direction. In a pertur- 
bative treatment the angle 9* is taken with respect to the direction opposite 
the direction of motion (Fig. 13). If only small angular components along 
the direction of the force are considered 

I{9*) (x cos 2 9*. (22) 

The intensity is essentially constant at small angles in the rest frame. 

The relation between the energies and angles in the lab and radiating 
particle rest frames is given by 

* - ?f, (23) 

« - M 

9 = (25) 

The direction of the radiation in the radiating particle rest frame is at CESR, 
for example, 9* ~ 0.03 or two degrees away from the backward axis. 

In the radiating particle rest frame the cutoff frequency is inversely pro- 
portional to the duration of the perturbation, which is cr*/c. Using equations 
24 and 25, and the relativistic formula length-dilation, the following relations 
are obtained 

a* z = a z h, (26) 
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< ~ 0(£) (27) 

^ ~ °(^2)> ( 28 ) 

which shows that the cutoff frequency at large angle does not depend on 7, 
the first prediction by the Coisson model. At CESR, lo c ~ 10 16 sec _1 , which 
is of order of the visible light frequency. 

The polarization vector of the emitted radiation in the radiating particle 
rest frame is given by Bj 

EW = i n,X(n ' Xn (29) 
where n* is the unit vector along the direction of observation. 

Using Equation 29, and the condition of orthogonality between E*, B* 
and n*, the three vectors are 

E* = K(6* 2 sin cos <p,0* 2 sin 2 <f>- 1,-6* sin <f>) (30) 
B* = K(-l + e* 2 /2,0,-e*cos(/)), (31) 
n* = (0*cos0,0*sin0,-l + 0* 2 /2), (32) 

with K a constant. The polarization component along x and y in the labo- 
ratory frame are 

E x = i(E*-B;) = + K^-sm2ct) = +K-^ sin 20 (33) 

E y = 7 (^ + j B*) = -kI^! C os20 = -^-^cos20. (34) 

Thus each component has four azimuthal zeros, and information is replicated 
every 45 degrees, which is the second prediction of the Coisson model. 

The total energy radiated in the laboratory frame can be expressed as 
an average over the boosted photon energies in the rest frame, times the 
number of photons N 

U = J2 l( k * + K) = N < ~/{k* + k* z ) >= N~f < k* > . (35) 

The energy flowing into a detector covering a solid angle cZfi, located at 
large angle, can be easily computed in the radiating particle frame. Using 
1(9*) ~ 1, Equations 23 and 25, and neglecting factors of order one an 
expression for the large angle spectrum is obtained, 

AU ~ N6*A9*A<j) J 2 = -^A#A0. (36) 
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The energy in the lab frame, U, contains a dependence on j 2 . The angular 
factor integrates to a constant (which agrees with Equation 8 in Ref. ||), 
leaving the I/7 2 dependence. This is purely due to kinematics. At CESR, 
for example, lOnW of visible beamstrahlung are available between 6 and 7 
mrad. 



Appendix B 

A beam-beam interaction simulation was developed from the program de- 



scribed in reference [10]. Gaussian beams in all three dimensions are as- 
sumed. Beams are sliced in 3-dimensional cells. The cells are typically 
0.25-0.5(7 along each axis and extend out to 3-4<r in each direction. Thus a 
total of 10 3 to 3 x 10 4 cells are simulated. The beams are then made to cross 
each other. In the first step, the first layer of the positron beam encounters 
the first layer of the electron beam. The electric fields are purely transverse 
to 0(1/7), and are computed assuming that the charge is located in a sphere 
located in the center of the cell. This is the "cloud-in-cell" model. Assuming 
cylindrical coordinates, a cell in beam one gets a total transverse deflection 



Ar lj = — (37) 

The summation runs over the cells in the opposite layer, b^- is the impact 
parameter between cell j in beam one and cell i in beam two, and P21 is 
the fraction of charge in cell i. At the end of each layer-layer interaction 
positions and velocities are updated, 

r'j = r'j + Ar'j, (38) 
Tj = Tj + r'jAz. (39) 

Az is the unit step taken along the beam direction. This allows for dynamic 
beams, with each beam pinching the other as the collision progresses, and 
the luminosity is computed as an overlap of the dynamic density functions. 

The program of reference [flOf] was found to be unfit for the simulation 
of flat beams. If the lattice is chosen to have the same number of cells in 
each dimension, the cells will be as flat as the beam. If the charge is then 
concentrated in the centers, a large force will be calculated, where in reality 
the total force is small, due to the cancellations of the large x components 
in the integral over the cells. Figure 14 illustrates this. 
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To reduce this problem the number of cells in x should be enlarged to 
make each cell square in the transverse plane. This solution is very CPU 
time-consuming. A solution was found by replacing each cell with a line of 
charge, called a "matchstick," and computing the integral 

Ar'ij = — ^2r e Pi I ^ x ^ x }^j_ ^q-j 
7 J b ij 

For the purpose of improving the convergence of the program, the match- 
sticks were kept horizontal throughout the interaction. Assuming matchstick 
lengths Li and Lj, the solution to the integral above is 

—2N2T P 4 4 
Ar'ij = ^-^(52 tnfn + 2b y g n , V 2t n g n - b y f n ), (41) 

where 

7" . _l_ T . 

(42) 
(43) 
(44) 

(45) 

(46) 
(47) 

Given the deflection vector, r', the total radiated energy with both x— and 
y— polarization is computed using 

F = ^Ar'. (48) 
2Az y ' 

The energy vector U for each beam is computed by summing 

2N 

U x = "£AU xj = £ — P^fFjAz, (49) 

2N 

U y = J2 AU yj = E^W^- (50) 

The program continues to interact the beams, layer by layer, updating 
trajectories with Equations 38-39, until the beams fully cross each other. An 



h 


= b x 


1 L i + Lj 
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= b x 
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Quantity 


Bin= 0.25o- 


Bin= 0.3a 


Fit 


Analytic 


Dyn. beams 


U x (W 12 eV) 


.4002 


.3979 


.4055 


.4051 


.4088 


U y (W 12 eV) 


.4013 


.3997 


.4049 


.4051 


.4163 


L/L 


1.00 


1.00 


1.00 


1.00 


1.12 



Table 2: Comparison between two different binnings, the fitted values, ac- 
cording to Equations 51-52, and the analytic predictions, all for stiff beams. 
The last column shows the same quantities for dynamic beams. 



option was inserted in the program to use or not to use Equations 38-39, that 
is to make the beams either dynamic or stiff. The reason for the option was 
to compare against existing analytic predictions for beamstrahlung given in 
reference ||. 

Appendix C 

The simulation program described above is used to make comparisons with 
the analytic predictions of reference || which are valid only for stiff beams. 
A slow, quadratic convergence was found when diminishing the cell size. 
When a cell of one beam overlaps with one from the other beam, the program 
computes a zero field. In reality, the contribution of nearby particles is 
important, due to the 1/6 dependence of the field. Because the emitted 
power depends on the field squared, the dominant convergence is quadratic. 

To adjust for this fact using finite computer resources two different lattice 
sizes a and b were used. The "exact" emitted energy U ex was extracted using 
the linear system 

U a = U ex -aa 2 , (51) 
U b = U ex -ab 2 , (52) 

and solved for U ex and a. 

The stiff-beam comparison of U, with different cell sizes and against the 
analytic predictions of reference || using the beam parameters of Table || 
are shown in Table |2[ There is agreement between analytic and simulation 
at the 0.2% level. If the beams are dynamic particles will move during the 
collision by about 1 micron, leaving the beam virtually unchanged in x but 
generating a substantial (O(10%)) squeezing in y. The squeezing will have 
two effects: it will increase the luminosity and it will generate slightly more 
power. There will also be a slight asymmetry between x and y and some net 
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polarization. The luminosity increases by 12%. The luminosity calculation 
was checked, for round beams, against the program of reference and 
our simulation agrees to within 1%. Figure 15 shows the analytical versus 
simulation comparison of U y /U x when two flat beams are separated by a 
vertical offset. We conclude that our simulation method has a precision of 
order few per thousand for beamstrahlung computations. 
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Figure 1: A general beam-beam collision. Seven parameters can be seen, 
corresponding to two transverse dimensions for each beam, a two dimen- 
sional impact parameter vector connecting the two beam centers, and one 
relative rotation in the transverse plane. 
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Figure 2: Normalized power emitted in beamstrahlung, as a function of 
normalized y— offset, a) e = 0.02. b) e = 0.04. The distance from minimum 
to maximum is shown, in units of a y . Uq is defined in Section 4. 
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Figure 3: The four beam-beam pathologies that lead to wasted luminosity; 
a) a y— offset; b) y— bloating; c) x— bloating; and d) a beam-beam rotation. 
The pathological beam is represented by the dashed ellipse. 




Figure 4: The beamstrahlung diagram corresponding to a perfect beam- 
beam collision. The two vectors are exactly equal. The dashed arrow is 
slightly displaced for display purposes. 
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Figure 5: Beamstrahlung diagrams corresponding to the four pathologies of 
Figure 3. The tips of vectors in part a are displaced for display purposes. 
Stiff beams are assumed. 
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Figure 6: Beamstrahlung diagrams for the same conditions as Figure 5, but 
assuming dynamic beams. 
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Figure 7: Beamstrahlung diagrams for the same conditions as Figure 5, but 
assuming an x— offset of 0.06a x . 
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Figure 8: Functional dependence of the beamstrahlung asymmetries defined 
in the text versus the waste parameter of Equation 3. 
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Figure 9: The six possible configurations arising from combinations of any 
two of the pathologies of Figure 3. a) y— offset and y— bloating, b) y— offset 
and x— bloating, c) y— offset and beam-beam rotation, d) y— bloating and 
beam-beam rotation, e) x— bloating and beam-beam rotation, f) y— bloating 
and x— bloating. 
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Figure 10: Beamstrahlung diagrams corresponding to Figure 
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Figure 11: Beamstrahlung diagrams, corresponding to Figures 9 and 10, 
after correction of the dominant asymmetry. Compare with Figure 3. 
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Figure 12: The dependence of the first asymmetry A\, as defined in the text, 
versus the vertical offset for the case of a vertical offset plus a rotation. 
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Figure 13: Dipole radiation in a radiating particle's rest frame. Indicated 
are the direction of the force and the angle corresponding to the observation 
angle in the laboratory frame. 
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Figure 14: Cell-cell interaction in the simulation program. The cell has an 
aspect ratio similar to the beam aspect ratio. In the "cloud-in-cell" model, 
all the charge is concentrated in a point in the center of the cell. In the 
"matchstick-in-cell" model, the charge is spread over a line along the cell. 
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Figure 15: Radiation polarization versus beam-beam offset. The solid line 
is the analytic prediction from reference ||, and the dots are from the sim- 
ulation described in the text. 



